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ABSTRACT 

We examine whether a cosmologically significant distribution of dark galaxy 
group or cluster-sized objects can have an optical depth for multiple imaging 
of distant background sources which is comparable to that from known galaxies 
while at the same time producing angular splittings of the same order of magni- 
tude. Our purpose is to explore whether such objects could realistically account 
for some of the observed lenses. Modeling such systems as isothermal spheres 
with core radii, and assuming a Schechter-type distribution function, we find that 
independent of the cosmology (open, flat matter dominated, or flat cosmologi- 
cal constant dominated) an allowed parameter range exists which is comparable 
in velocity dispersion to that for known compact groups of galaxies, although 
the preferred core radii are somewhat smaller than that normally assumed for 
compact groups. Dark cluster-sized objects, on the other hand, cannot repro- 
duce the observed lensing characteristics. If the one known Dark cluster were 
a good representative of such a distribution, most such objects would not pro- 
duce multiple images. We also present a result for the angular splitting due to 
an isothermal sphere lens with non-zero core radius, extending earlier work of 
Hinshaw and Krauss (1987). Our results are expressed as contour plots for fixed 
lensing probabilities, and angular splittings. 
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1. INTRODUCTION 



Ever since the first multiply imaged quasar was observed, it was clear that the statistics 



of gravitational lensing could be utilized to pin down cosmological parameters (i.e. Turner 



Ustriker, fc Gott 1984| ; |Kundic et al. 199^j ). Specifically, for a given cosmological model, one 



can predict the optical depth due to lensing by normal galaxies (presuming one has a model 
of normal galaxies and their distribution), and compare that with observations. Several 
large scale surveys have been performed, searching for multiple imaging of distant quasars 
by intervening galaxies, and more than a dozen such lensing events have been observed 
( [Kochanek 1996|) . However, there is one slight peculiarity. While the overall frequency 
of lensing events, and the rough angular splittings are reasonably consistent with model 
expectations, in a significant fraction of the cases, the actual lensing galaxy is not visible. 
Given the predicted mass of the lensing systems, it is not obvious that such a large fraction 
should remain unresolved. 

This prompts the natural question: Could the lenses be dark objects, perhaps objects 
related to the distribution of dark matter in the Universe — perhaps failed galaxies? This is 
not a new idea, and it is one which has been beset with problems. In the first place, anything 
close to a closure density of compact objects generally produces an optical depth for lensing 
which is too large QPress fc Gunn 1973|) . Second, if the dark objects are failed galaxies, then 
their properties will generally preclude producing multiple images with the observed angular 
splittings ( |Hinshaw fc Krauss 1987| ). 



Here we examine another possibility. Could larger systems, on the scale of groups or 
clusters account for an observable fraction of the known lenses? The recent observation 
of a cluster-sized mass distribution containing one luminous galaxy flHattori et al. 1997| ; 
[Mushotzky 19971 ) provides some additional a posteriori motivation for considering this hy- 
pothesis. 

On first glance such a possibility seems implausible. Larger systems, with larger veloc- 
ity dispersions, will produce larger angular splittings, if they produce multiple images at all. 
Thus, it seems unlikely that such systems might reproduce the observed lensing character- 
istics, which, as alluded to above, are comparable to those one would predict for the known 
population of galaxies. However, Hinshaw and Krauss (1987; hereafter HK) demonstrated 
that under certain conditions a finite core radius suppresses the mean angular splitting due 
to isothermal sphere mass distributions. Here, we generalize the earlier HK result, and prove 
that this approximation is good for all lensing impact parameters inside the critical disk for 
multiple image formation. We then demonstrate that for a reasonable range of velocity dis- 
persions, core radii, and total mass, assuming a Schechter-type distribution function, such 
systems can produce optical depths for lensing comparable to that due to known galaxies, 
with comparable angular splittings. 



It is well known that due to a variety of selection effects (magnification biasing, etc.) 
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the actual fraction of strong lensing events in any sample can differ dramatically from the 
naive optical depth calculation. However, because these systematic effects should be largely 
independent of the nature of the lenses themselves, if they are producing comparable angular 
splittings, etc, we need not consider these effects in detail here. In particular, we compare 
the calculated optical depth for lensing by dark objects to the naive optical depth for lensing 
by known galaxy distributions, known to be in the range of 10~ 3 to 10~ 4 . Presumably, if the 
other selection effects are comparable, this will then imply that such systems could produce 
at least some fraction of the observed events. We thus derive contour plots in the parameter 
space of velocity dispersion and redshift for fixed total optical depth in order to explore 
the suggested range of dark lenses. For this range we then explore the magnitude of the 
predicted mean angular splittings. We conclude with a brief summary of our results. 



2. THEORETICAL FORMALISM 



2.1. Optical Depth for Lensing 



We start with a mass density distribution with the following form: 



p(r) 



1) 



2nG(r 2 + r c 2 ) 

where a is the velocity dispersion of this system, and r c is the core radius. In order to derive 
a finite total mass for the system we must assume some cutoff radius r M . The total mass is 
then given by: 



m(< r^) 



2fL 
~G 



roo - r c tan 1 ( — 



(2) 



If we assume a matter density fraction Q Q in the Universe, then we can calculate the 
number density of objects, n, as this fraction of the critical density p crit = divided by 
the mass of each system: 



n = Q Q - 



Pcrit 



3H G 2 



"oo - r c tan 1 (^) 



Following Turner, Ostriker, & Gott (1984) we then define a lensing probability factor 



(3) 



c 3 imal 
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(4) 



Here, ao = 47r(-) 2 , is the bend angle for an isothermal sphere with core radius r c = 0. 



Then, the differential optical depth is: 
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where y oi is the angular size distance ( |Peebles 1993| ): 



H o a ri 1 / r— /■* dz , 

Voi = = smli U/S2r / (6) 



and 

/(/%) = 1 + 5/3 - -/3 2 - - v^(A) + 4) 1 (7) 

We will explain the origin of the function /(/3o) and the parameter /5q in the next subsection. 
In the above equations, ao is the scale factor of the Universe, r, is the coordinate distance, 
Zi is the redshift, y g is the distance between the observer and the lensed galaxy, y^ s is 
the distance between the lens and the source, and y os is the distance between the observer 
and the source. f^R is the curvature term, and is the cosmological constant, so that 
Q Q + Q-r, + Q\ = 1 in all cases. 

In a flat universe, ye s is simply y os — y oi . However, in an open universe (with = 0): 



Via = Vos V 1 + ^RVol ~ Voi A/ 1 + ^RVos 



If all lensing systems had the same mass and core radii, this would then be sufficient 
to calculate the relevant optical depths. However, if one has a distribution of masses and 
core radii, it is necessary to integrate over this distribution. For simplicity, we assume that 
these systems are distributed, like observed galaxies and clusters, with an effective Schechter 
function, with a parameter L, which is related to the total mass of the system. 

mdL=(±Yew(-±)d(£) (9) 



L*J V L*J \L 

We can then consider a and r c as functions of L, and then integrate over this distribution 
when deriving optical depths. We use the conventional parameter L to parametrize this 
distribution. However, because we are interested in dark structures, L should be understood 
to refer to mass, rather than luminosity in this case, unlike its representation for luminous 
objects. Although one could easily choose to other arbitrary mass distribution functions for 
equation (|]), it seems reasonable to use the Schechter function, since this fits the distribution 
of luminous galaxies. Our chief purpose here, in any case, is to explore what general mass 
ranges are picked out in order to fit the data, and this rather general parametrization is as 
useful as any other for this purpose. 



2.2. Lensing Cross Section and Angular Splitting 

In order to determine the cross section for strong lensing events, one needs to determine 
the bend angle, a, of the light trace from the source, which in the case of an isothermal sphere 
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with a finite core radius, is a function of ( Hinshaw fc Krauss 1987 ) the velocity dispersion 
a, impact parameter b, and core radius r c . The general formula for the bend angle is: 



a(b) 



46 r°° , <9$ 
dr- 



c 2 Jb 



where 



Or i/r 2 — b 2 
$(f) = - / d 6 ? 



]3j G P( r ') 



(10) 



(11) 



HK derived a closed form approximation for a using a mass distribution of the form. 
Here, because there will be a relation between and the number and mass density of dark 
objects, we have to refine this earlier result to allow for finite r M . It is simple to show that 
in this case, a first order expansion in {r 2 + b 2 )/r' 2 <J leads to 



b 2 + r 2 — r c 



a(b) = ao 



Following HK, we define A = a CT /ao, where 



'b 2 + ri 



7T7Y 



(12) 



so that 



c a yoi yes 
H G (1 + zi)y 

b + l 



A 



a(b) 



(13) 



(14) 



where / is the transverse distance of the lens center from the line of sight. 
Defining 

1 + — 



7T7V 



one can show 



r r = 



b + I = a c 



1 + 



7T7V 



b 2 + r 2 — r c 



(15) 
(16) 

(17) 



This equation is identical in form to that derived by HK, so we can then use their results to 
directly write down the cross section 



a cs = 7ia 2 c J((3 ) 



where 



A) 



(X rr 



< - 



(19) 



Also note that /(/?o) is zero when (5q is larger than 1. 
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The bend angle not only allows us to derive the lensing cross section, but also the 
observed angular splitting between lensed images. HK demonstrated that for a lens located 
along the line of sight to the source, this splitting is given by 2^/1 — 2/3 ctQ— . In an 
appendix we present an analytic proof that this approximation is good to 10% for all values 
of I relevant for multiple image formation. Hence we include this factor in our analysis. 
Note that, depending on the value of fl , this can lead to significant reductions in the mean 
predicted splitting, allowing systems with velocity dispersions larger than those of galaxies 
to produce comparable angular splittings. 



3. NUMERICAL RESULTS 

Utilizing the formalism developed in the last section, we can estimate the optical depth 
due to dark objects of specified velocity dispersion and core radius for several different 
cosmologies. We consider first the results obtained without averaging over a distribution of 
objects, and then explore how such averaging can impact upon the optical depths obtained. 
We then focus on considerations of the mean angular splitting produced by such dark lenses, 
and finally make some comparisons to observations. 

Figure 1 displays contour plots for optical depths of 10~ 3 (thick lines) and 10~ 4 (thin 
lines) with no averaging over a distribution of masses, but instead taking a uniform cosmic 
density of objects of fixed mass, leading to net density Q - Note that there are three free 
parameters, a, r c and r M . For a fixed ratio y 2 -, related to the mass per object, the integrated 
optical depth out to source redshift z s is then simply a function of c 2 /r c , so it seems reason- 
able to present the results in the phase space for this ratio, normalized in units of ^prff^ 
vs z s . Comparison to existing, or proposed, mass distributions can be done by plugging in 
physical values for a, r c , or both. Finally, all the black curves correspond to y 2 - = 10, and all 
gray curves correspond to ^ = 100. The solid curve sets are for a flat universe, with Q = 1, 
the dash-dash curve sets are for a flat universe model with f2 = 0.3, and = 0.7, and 
the dash-dot curve sets are for an open universe model with Q Q = 0.3 and zero cosmological 
constant. 

Note first that the — = 100 curve is higher than the — = 10 curve for a fixed cosmo- 
logical model. This can be understood as follows. For a fixed value of r c and a, increasing ^ 
increases the mass per system. However, since the overall number density of systems is de- 
termined by the requirement of some fixed fl Q , this number density then decreases inversely 
with leading to fewer lenses. Since the lensing probability for isothermal spheres (with 
core radius) is not determined by their total mass, but rather by their velocity dispersion 
(and r c ), having fewer lenses means a smaller optical depth, even though the mass per lens 
increases. Next note that while it is known that for a fixed density of lenses, the optical 
depth increases for an open universe model compared to a flat one fl — l, this effect is not 
seen here simply because we normalize the number density by the matter density. Since an 
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open universe has a lower matter density than a flat matter dominated universe, the number 
density of lenses decreases in proportion, and this is reflected in the optical depth contour 
plots shown here. The well known fact that the optical depth for a cosmological constant 
dominated flat universe is larger than for a matter dominated flat universe is also reflected 
in our results. 

Before proceeding to compare with observational data to determine if these estimates 
might correspond to reasonable structures, we need to allow for the fact that the dark lenses 
are distributed over different masses as per equation (|9|). 

Without strong theoretical guidance in this regard, we choose, by analogy with luminous 
objects, «i = —1 in equation @. We also assume % is proportional to L 012 , where we choose 
a?2 to be negative to avoid a singularity. We display our results for ct2 = —0.5 and —1.5. 
Note that in their lensing probability analysis for luminous objects, Krauss and White (1992) 
suggested that a<i could be a small positive number, although not much larger than zero at 
best. Note also, however that if one assumes the core radius is not a function of mass, then 
a?2 is negative. In Figs. 2 and 3 our results for the two different choices of a<i are shown. As 
can be seen, the curves move toward smaller values of ^*/{h-^Mpc) wnen a 2 decreases. That 
is, the probability is higher for a given source redshift and ^/{h-^Mpc) > wnen a 2 decreases, 
as one might expect, since this implies that larger core radii are less heavily weighted. We 
should remind our readers here that a* and r* are normalized to L* in equation ([|). 

We next consider the mean angular splitting induced by the lensing distribution. Again, 
we first consider the case when we do not integrate over a mass distribution of lensing. 
Recall that for a singular isothermal sphere distribution, the mean angular splitting is a®. 
In the case of finite core radius, the expectation value of — , using our result that A9 pa 
2VT=W a f-, is: 

yos 

<x>= — J dr x 

A9 1 f 1 A9 

< >= - / dr (20) 

ao r J ao 

where r is the optical depth from equation (||). In Fig. 4, we show contour plots for the case 
when the expectation values of the angular splitting either 0.5 or 0.8 times that expected 
for the equivalent singular isothermal sphere distribution. The thicker curves represent the 
expectation value of 0.8a , and the thinner curves represent the expectation values of 0.5a . 
Other symbols are the same as in Fig. 1. The vertical axis is still ^/(h-^Mpc) an< ^ ^ ne horizontal 
axis is source redshift. 

Note that the ^ = 10 curve is higher than the ^ = 100 curve. This is opposite to 
the probability plots discussed in the earlier subsections, and reflects the simple fact that 
larger r c for a fixed a results in a smaller angular splitting. This tension between angular 
splitting and optical depth is important, because it implies that to keep angular splittings 
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small enough to be comparable to those observed, one cannot allow too large a value of a 
and still keep the overall optical depth comparable to that due to luminous galaxies. This 
tuning is an important feature which constrains this scenario. 

Finally, we repeat the above analysis, including an integration over a Schechter distri- 
bution of lenses, as described above. The parameters are the same as discussed above. In 
Fig. 5, we chooose = —0.5, and the thick and thin curves represent when the expectation 
values of the angular splitting are respectively 0.9 and 0.7 of the zero core radius case. In 
Fig. 6, we choose «2 = —1.5 and the thick and thin curves represent when the expectation 
values of the angular splitting are 0.85 and 0.95 of the zero core radius case. Note that 
when a?2 decreases, the expectation values for the angular splitting approach the isothermal 
sphere result, as one might expect, since for these values, the core radius is becoming less 
significant. Thus, if a 2 is less negative than we have assumed here, the expectation value 
of the angular splitting could even be smaller in comparison to the equivalent isothermal 
sphere distribution than we have found. 



4. COMPARISON WITH OBSERVATIONS AND CONCLUSIONS 

Over a dozen multiply imaged quasars have now been observed in various optical and 
radio surveys. While these have angular splittings characteristic of those one might expect to 
be induced by galaxies (i.e. 0(2" -?>")), several have no observable lensing systems, even when 
one might expect the lensing galaxy to be resolvable. When one calculates the optical depth 
for lensing by known galaxy distributions, ignoring magnification biasing, one typically finds 
10~ 4 < t < 10~ 3 . Thus, one might expect that a predicted optical depth of as low as 10~ 5 for 
some other distribution would produce at least one or two lenses in the existing surveys. In 
fitting the optical depth, however, one must confront the tension between optical depth and 
angular splitting discussed above (also compare Fig. 1 and Fig. 4). Systems with velocity 
dispersion greater than about 450km / sec would require the expectation value of the angular 
splitting to be less than half that predicted for an equivalent singular isothermal sphere. 
This, in turn, requires r c to be larger relative to a, which in turn, however, suppresses the 
optical depth. Nevertheless, in Table 1 we display several sample values of core radii and 
velocity dispersions which would be expected to produce a fraction of the observed lenses 
with angular splittings comparable to those observed. In this table we have fixed ^ = 10, 
set the source redshift to 3, and fixed the parameters so that the mean predicted angular 
splitting is 2.5 arcsec. Note that the parameter range is comparable to compact groups of 
galaxies, which have compatible ^/(i^^Mpc) • ^ OT exam pl e ; °ne study suggests a is around 331 
km/ sec, and r c is around 15h~ 1 kpc ( |Mendes de Oliveira fc Giraud 1994j ; [Pildis, Bregman] 



|fc Evrard 1995| ). These numbers give ^( CT A"V^ C ) ~ 7 3 * iq 6 , as compared to the preferred 



values of j?j(™i S Mp C ) between 10 7 and 10 9 to result in an optical depth comparable to that 
of the known galaxy distribution. Also note that if we were to integrate over a distribution 
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of such objects, our earlier arguments suggest that comparable optical depths and splittings 
could be obtained even if the mean value of r c were somewhat larger than given in the table. 

However cluster-sized objects do not seem to be viable candidates for dark lenses be- 
cause of the larger predicted splittings when multiple image formation does occur, and more 
importantly because the predicted optical depth is too small. (The r c of a regular cluster is 
large, but the velocity dispersion is only about 2.5 times larger than that of compact groups 
of galaxies flZabhidoff et al. 1993| ). ) This conclusion is reinforced by the recent observation 
of an actual dark cluster by Hattori et. al. (1997) in the lensing system, MG2016+112. 
Based on the inferred mass of the object using X-Ray estimates of the potential, its size, and 
the size of the core radius of this system, we find that (3 > 1/2, implying that the cluster is 
not responsible for the observed multiple images in this system. This is again supported by 
the small angular splitting of 3.4 arcsec between the images, which is characteristic of the 
one observed galaxy in this system. 

Thus, perhaps paradoxically, the only known example of a large-scale dark object sug- 
gests that such cluster-scale objects, even if they have a significant mass density in the uni- 
verse, are probably largely irrelevant for the statistics of multiply imaging distant quasars. 
Rather, dark objects on intermediate scales, between galaxy and cluster scales, are more 
likely possible candidates for dark lenses. A very recent paper which has shown a very high 
rate of galaxy lensing in radio surveys may be significant in this regard flJackson et al. 19981 ) • 
We also note both that a recent study of the luminosity function of the compact groups of 



galaxies also lends support to our choice of Schechter ot\ = — 1 QZepf, de Carvalho, fc Ribeiro 



1997|) , and that the currently favored cosmological model involving a flat universe with cos- 



mological constant (i.e. [Krauss 1998| ) produces the largest optical depths at high redshift for 
given ~/7pr|^L ) as one would expect, based on the increase in the optical depth for lensing 
by galaxies in this cosmology. 

Clearly, in order to know whether dark clusters are important for lensing in the actual 
universe, larger surveys will be required, in order to reliably determine how many lensing 
events might not be associated with galaxies. If such events continue to be observed, our 
analysis suggests that the distribution of angular splittings will be an important observable 
which might constrain possible models. 
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A. APPENDIX 

We prove here the angular splitting, A9, is a monotonically decreased function of L = 
(not to be confused with the quantity L in the Schechter distribution function used in the 
text), where / is the distance from the line of sight of the lens. Then, using this fact, we 
can show that A9 is well approximated by 2^/1 — 2/3 «o^, the angular splitting derived by 
Hinshaw and Krauss for lenses along the line of sight. We start by redefining the parameters 
of equation (14) in Hinshaw & Krauss (1987): 



x = — (3= — 



(Z cr 



(Al) 



We then find 



x 3 + 2Lx 2 + (L 2 + 2f3- l)x + 2Lf3 = 



(A2) 



(We remind the reader that (3 is between and |, and \L\ < 1.) The above algebraic equation 



is solvable. Setting x = y — ^ the solution is 



y = U + V, ~(u + v)+i^Y-(u-v), 



■-{u + v) -i— (u 



(A3) 



with 



u 



\ 



\ 



Q 
2 



(A4) 



-p = 1 - 2(3 + -L 2 > 
3 

1 L 

2 3 



l + (3-lL 2 )>0 



(A5) 
(A6) 



If we want to have three real solutions of y, which correspond to three lensing images, then 
u — v should be a pure imaginary number. However, u + v is a real number, and this gives 



us u* = v, and J + (^j is a pure imaginary number. We thus define 



I 11 



and 



r 2 = ±-[(1 - 2[3) 3 - L 2 {2 + 10/5 - (3 2 ) + L A ] > 

— f 



r = 0^L 2 = ^ = f{(3) 



(A7) 
(A8) 

(A9) 



so we can re-write u and v as following: 



u = 



6 q 



+ r 2 exp(— i(f>/3), 



v = 



+ r 2 exp(i0/3) 



(A10) 
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where 



this implies that 



Recall that r 2 + (f ) 2 



2r 

tan0 = — > 
Q 



< < 



7T 



-p\3 



, then we have 



u — — \ 



'-p 



exp(-«0/3), 



v = -\ 



-P 



exp(i0/3) 



(AH) 
(A12) 

(A13) 



we can write down the three solutions of y: 



u + v 



-\{u + v) +i^(u-v) 



-2^cos(f) < 
2 v /^cos(f + f) > 

2 v4 cos (f-f)>^ 



yi = 

V2 = 

= -\{u + v) -i&(u-v) 
therefore, the angular splitting is: 

^ = a (A(x + L)) = a Ay = a (y 3 - y±) = 2a y/^pcos - ^ J > 

y os \ / 

Based on the above equation, one can derive (after some work) that: 
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dL V~ V ~dL, 



< 



(A16) 



Note that AO varies from L = to L = L 

I J 

y os 



AO = 2Jl-2(3a 



when L = 



AO = V3\ 1 - 2(3 + -Ll a — , when L = L c 
V 3 y os 



(A17) 



Finally, the average ratio of AO to 2^1 — 2(3 a ^- at L = L is 



r dp 

Jo 



AO 



0.88 



(A18) 



Since we have shown that AO is a monotonically decreasing function of L, this implies 
that over the relevant range of L, AO does not differ by more than 10% from its value at 
L = 0. 
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TABLE 1 

Sample Values of Core Radii and Velocity Dispersion Relevant to Lensing 





model 




a (km/ 'sec) 


r c (h 1 kpc) 


r 


<A6> 




= i,n A 


= 


427 


4.6 


2.9* 10" 5 


0.47a 


n = 


0.3, n A 


= 0.7 


369 


5.0 


6.9* 10" 5 


0.55a 


n G = 


0.3, n R 


= 0.7 


416 


4.4 


3.4* 10" 5 


0.50a 



Note. — The values are obtained by assuming ^ = 10, source redshift at 3, and < A6 > is fixed at 2". 5. 
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Fig. 1. — Contour plots of lensing probability for fixed lens mass. The thicker width lines 
represent an optical depth of 10~ 3 , and the thinner width lines represent an optical depth of 
10~ 4 . The black curves correspond to ^ = 10, and the gray curves correspond to ^ — 100. 
The solid curve sets are for a flat universe model, with Q a — 1. The dash-dash curve sets 
are also for a flat universe model, but with Q = 0.3, and = 0.7. The dash-dot curve sets 
are for an open universe model, with Q Q = 0.3 and zero cosmological constant. 
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Fig. 2. — Contour plots of lensing probability as for Fig. 1, but integrating over a Schechter 
distribution of lenses 
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Fig. 3. — As for Fig. 2, but with a 2 = —1.5. Note that all curves involve smaller value of 
f'/jHtfi) compared to those in Fig. 2. 
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Fig. 4. — Contour plots of the expectation values of angular splitting normalized to that 
for an isothermal sphere, for fixed lens mass and showing contours < ^ >= 0.5 and 0.8. 
The thicker width lines represent < >= 0.8, and the thinner width lines represent 
< — >= 0.5. All other features are as for Fie;. 1. 
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Fig. 5. — As for Fig. 4, but integrating over a Schechter distribution, with = —0.5 and 
showing, < ^ >= 0.7 and 0.9. 
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Fig. 6. — As for Fig. 5, but for ct2 = —1.5 and showing, < ^ >= 0.85 and 0.95. Note that 
all curves involve smaller values of ^,!n^itT\ compared to those in Fig. 5. 



